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TMA MST204 03 Cut-off date 7 May 1998 


This assignment covers Units 9, 12, 14 and 15. 


Question 1 (Unit 9) 
(i) By using the Gaussian elimination method, determine whether the following 
sets of simultaneous equations have 
¢ no solution, 
e a unique solution, 
¢ an infinite number of solutions. 
In each case, if the set of equations has a unique solution find the solution, 
and if it has an infinite number of solutions find the general solution. 
(a) 32, + 522 — 6x3 4 
@ + 2tg + 32g 
3a, + Gro + Bag 
(b) 8a, + 5a — 6x3 
a + 2g + 3xy 
Bay + 6x + Ixy 


nueu 


“oui 
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(ii) To investigate the absolute conditioning of the set of simultaneous equations 
in part (i)(a), the value of the right-hand side in the third equation is changed 
from 8 to 8+ ¢, where € is a small number, so that the equations become 


3x, + 52 — 62g = 4, 
a + 2% + 3xy = 3, 
3a, + Gry + Bry = BtE, 


Find the exact solution of the modified simultaneous equations, What conclu- 
sions (if any) can you draw about the absolute conditioning of the original set 
of equations from these results? 


Question 2 (Unit 12) 


This question is concerned with steady-state heat loss through an external wall of a 
room of a house. The wall is a cavity wall, consisting of two layers of brick with an 
air gap between them. The wall is 2.5 metres high and 3 metres wide. The layers 
of brick are each 0.11 metres thick, and the width of the air gap separating them is 
0.05 metres. In the wall there is a window 1.5 metres high and 1.5 metres wide; it 
consists of a single pane of glass 0.006 metres thick. The thermal conductivity of 
brick is 0.6W m7! °C", and the thermal conductivity of glass is 1.0 Wm! °C>!. 
The convective heat transfer coefficients for the surfaces of both the brick wall and 
the window are 10W_m~?°C~! for the inside of the house and 20W m=? °C~" for 
the outside of the house. The combined heat transfer coefficient for the air gap 
of the cavity wall is 1.5 Wm~*°C~!. The outside and inside air temperatures are 
constant at 2°C and 20°C, respectively. 


(i) Calculate the U value for the cavity wall, and the U value for the window. 


(ii) Find the total amount of heat lost in one day through the cavity wall and 
window. 


Find the temperatures on each of the four surfaces of the brick layers of the 
cavity wall. 


(iii 


Suppose that the wall is to be insulated by filling the cavity with plastic foam, 
whose thermal conductivity is 0.034 Wm7!°C~!. What will be the resulting 
reduction in heat lost per day, assuming that all other conditions remain the 
same? What will be the temperatures on the surfaces of the wall inside the 
room and outside the house with cavity insulation in place, all other conditions 
being unchanged? 


(iv 


(U1) 


Question 3 (Unit 14) 


(i) 


(ii) 


Show that the area of a triangle ABC is equal to half of the magnitude of the 
cross product AB x AC. 


The position vectors a, b, of three particles A, B, C with respect to a fixed 
origin O are given by 


a=ti+j+(t+2)k, 
b = (2t+ 1i+ (t+ 1)j+ (t+ Ik, 
c= (t+2)i+ (t+4)j+ (2t+3)k, 


where ¢ represents time, and i, j, k are fixed Cartesian unit vectors, 


(a) Find the vectors AB and AC in terms of t, i, j and k. 


(b) Find a unit vector perpendicular to the plane containing the triangle 
ABC, and deduce that the particles move in such a way that successive 
positions of this plane are always parallel to each other, 


(c) Find the area of the triangle ABC. 


(d) Find the value of ¢ for which the area of the triangle ABC is a minimum, 
and find this minimum area. 


(e) Find the angle between AB and AC when the area of the triangle ABC 
is a minimum. 


Question 4 (Unit 15) 


a 


fi 45 
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A netball player throws the ball from a point h metres below and 44 metres hor- 
izontally from the net, with the purpose of making the ball fall through the net. 
The ball is thrown with a speed of 9ms~! at an angle of projection of a to the hor- 
izontal, Throughout the question you should neglect air resistance and model the 
ball as a particle. You should further assume that the acceleration due to gravity 
is constant and of magnitude g = 10ms~?. 


(i) 


(ii) 


(iii) 


(iv) 


From Newton's second law of motion, show that during its flight the horizon- 
tal and upward vertical displacements x and y of the ball from its point of 
projection satisfy the equation 

y= rtana— £1? sec* a. 
If the player throws the ball from a point = 1 metre below the net, find the 
possible angles of projection of the ball in order that the net should be on the 
ball’s trajectory. 


By examining the vertical component of the ball’s velocity, show that the ball 
actually falls through the net from above (rather than passing through the net 
from below travelling upwards) for both the possible angles of projection. 


Show that, for the given speed of projection and horizontal distance from the 
net, the maximum achievable height of the net above the point of projection 
for different angles of projection a is 2.8 metres. 


[You may find the following identity useful: 


sec” a = 1+ tan? a] 


[3] 


TMA MST204 04 Cut-off date 28 May 1998 


This assignment covers Unit M. 


Ther 


‘e two questions in this assignment, each marked out of 50 


The purpose of this TMA is to help you practise the modelling skills described in 
Unit M, which you will need when you come to tackle the modelling exercise for 
TMA 07. It consists of two questions. In the first we invite you to analyse a piece 
of modelling that has been carried out by a professional mathematician. We hope 
that this will help you to understand better the processes involved in modelling. 
This question is similar to those we asked in Section 2 of Unit M, about the drug 
therapy and the skid marks models. The second question is a problem for you 
to solve by mathematical modelling. and has a similar purpose to the exercise in 
Section 5 of Unit M. In both questions the organization of the work follows the 
outline of the modelling process given in the summary of Section 4. You may find 
it useful to glance over this to refresh your memory before you begin. 


This TMA is designed to take no longer than any of TMAs 01 to 03, Tf you find 
that you are spending a lot of time on it, then you are probably reading more into 
the questions than was intended. 


Question 1 


The following passage is part of an article describing a mathematical model of the 
costs and benefits of insulating houses against heat loss. Read the passage, and 
then tackle the questions which follow it, You will find that there are several gaps 
in the passage; some of the questions ask you to work out what is missing. However, 
you should be able to follow the sense of the article despite these omissions, 


Saving Heating Costs with Loft Insulation 


This article is concerned with the financial benefits which may be derived from 
insulating the loft, or roof space (the space between the ceilings of the uppermost 
rooms and the tiled or slated roof itself), of one’s ho This is done by laying 
insulating material such as fibreglass over what one might call the floor of the loft, 
so as to decrease the flow of heat through the ceilings of the rooms below. 


Loft insulation reduces heat loss from a house and saves the householder money by 
cutting the cost of heating; but against this saving must be set the cost of installing 
the insulation itself. Many houses in the UK lack loft insulation; and of those houses 
that have it, many have too little, The recommended thicknes for fibreglass loft 
insulation in the UK is 100 mm; few houses are insulated to this level. The cost of 
insulating one's house clearly increases if one installs thicker insulation, and this 
is no doubt a major factor influencing how much insulation people install; but 
using a thin layer of insulation may be a false economy if the extra savings on the 
heating bill from thicker insulation outweigh the additional costs of installation. 
The purpose of this article is to investigate how the overall ancial benefits of 
installing loft insulation vary with the thickness of the insulation installed, and to 
find the thickness for whicl’ the financial benefit is greatest. 


The cost of heating is a continuing, or recurrent, cost, whereas expenditure on insu- 
lation is a one-off, capital, cost. It is necessary to find a suitable way of combining 
these two rather different quantities. The method of combining a recurrent and 
a capital cost, which will be used here is to add to the recurrent cost the interest 
which would have been received if the capital had been invested instead of spent. 
The idea is that every £1 spent on insulation, if it had been invested instead at 
(say) 5% per annum, would have earned the investor 5p each year; by choosing to 
spend that £1 on insulation, one foregoes the 5p interest, and this can be regarded 


as an annual cost which accounts for the outlay on the insulation. 


The problem, therefore, is to find the optimum thickness of loft insulation, taking 
into account both the cost of heating and expenditure on the insulation. The 


approach to solving this problem which will be adopted in this article is to obtain 
an explicit formula for ¢(t), the overall annual cost in £ associated with insulating 
each square metre of the loft space with insulation t metres thick, and to find the 
value of t for which c(t) is a minimum, using the usual methods of the calculus, One 
consequence of this approach is that the only factors which are taken into account 
are financial ones, 


As has been explained above, the total cost ¢(t) is the sum of two components, the 
cost of heating, and the interest foregone on the cost of installation; each of these 
terms will also depend on the thickness of the insulation. Thus 


e(t) = h(t) + a(t), Q) 


where h(t) is the annual cost in £ of supplying the heat lost into the roof space 
through one square metre of the ceiling, and i(t) is the annual interest in £ received 
on an invested sum equal to the capital cost of installing one square metre of 
insulation, when the thickness of the insulation is t metres. Each of these two 
components will be considered separately. 


To find an explicit expression for h(t), it will be necessary to consider the loss of 
heat into the roof space of the house through the ceiling of each room immediately 
below it. In any particular case the rate of loss of heat will depend on the materials 
from which the house is constructed; but the model to be developed here will apply 
to any house whose ceilings adjacent to the roof are made from a single material 
and ave of uniform thickness, For simplicity’s sake, the effects of the timbers in 
the roof to which the ceilings are attached will be ignored; and only conduction 
(and not convection or radiation) will be considered. The temperature in the roof 
space will be taken to be equal to the outside temperature. The rate of loss of 
heat, q¢ watts, through 1 square metre of ceiling, when insulation t metres thick is 
installed, is given by 


q¢ =U (Bin — %) (2) 
= 5 (60 — Bout) (3) 


where U is the U’ value for the ceiling material alone, « is the thermal conductivity 
of the insulating material, and 6),, @ and @., ave the temperatures inside the 
house, at the interface between the ceiling and the insulation, and outside the 
house, respectively. Thus 


«U (Bin = Pout.) 
K+Ut ~ 
The heat that is being lost at this rate has to be supplied by the heating system; 
so the contribution of the losses into the roof space to the annual cost of heating 
the house is obtained by summing up the heat lost at this rate for the whole period 
of time that the heating system is operating. It is natural to assume that the 
required internal temperature 4, does not change through the winter months while 
the heating is on. The external temperature will vary, of course; but this can be 
allowed for by taking for oy. the average external temperature over the period in 
which the heating is on, Then the total heat lost into the roof space through each 
square metre of the ceiling during one year, Q Wh, is given by 
_ «UTAO 5) 
~ K+Ut" . 
where T is the total time for which the heating is running and A@ is the average 
temperature difference between the inside and the outside of the house. (The units 
for the total heat lost are watt hours. One watt hour, or Wh, is the energy produced 
by a source rated at 1 watt running for one hour. A more common unit of energy 
is the kilowatt hour, or kWh; one kilowatt hour is 1000 watt hours.) 


q= (4) 


Finally, the annual cost of providing the heat which is lost into the roof space 
through each square metre of the ceiling is obtained by multiplying Q by the cost 
of a watt hour of energy, say £C, which gives 


_ CxUTAO 


NST rare (6) 


The other component of the total cost is the cost of installation. The most eco- 
nomical way of insulating one’s loft is to do the job oneself; it is not a very difficult 
job, though a dirty and tiring one, so it seems not unreasonable to suppose that 
this is how it is done. It is then necessary to take into account only the cost of the 
insulating material, This is easy to calculate because it will be proportional to the 
volume of material required. Suppose that the price of 1 cubic metre of insulating 
material is £P; then 


IPt eg 
a(t) = 00" (7) 
where the rate of interest is 1%. Combining the two costs gives 
CrUTAO , IPt 
e(t) = h(t) + a(t) = Sane, + oo (8) 


It remains to find the value of t (if there is one) for which ¢(t) is a minimum, Now 


« = | deliberate omission! (9) 


thus in order for ¢(t) to have a minimum, t must satisfy 


(x + Ut)? = [deliberate omission! |, (10) 


and therefore 


/100CKTAO Kk 
t= cay Ieee ie (11) 


This gives the thickness of insulation for which the total annual cost is a minimum. 


As an example, consider the following data. 

¢ The ceilings are constructed from plasterboard which is 10mm thick and has 
a thermal conductivity of 0.16 Wm! °C>!. 

« The insulating material to be used is fibreglass; the fact that it is usually avail- 
able only in certain standard thicknesses may be ignored. A voll of fibreglass 
insulation 8m long, 0.35m wide and 100mm thick costs £4.80. The thermal 
conductivity of fibreglass is 0,04 Wim! °Co!, 

e Interest on invested money is 5% per annum. 

¢ Heating is required in the average house in an average year for 15 hours during, 
the day, for 180 days during the year, and the average temperature difference 
between the inside and the outside of the house in that period is 8°C. 


« The running costs of the heating system are 2.5 pence per kilowatt hour. 


The U value for the ceiling material alone is 16 Wm~?°C-!. The cost. of 1m® of 
fibreglass is £17.14. The value of t at the minimum is therefore 


TOO X00 XOUEXISOXISXS 0.04 
V 1000 x 5x 17.14 iris 


and the optimal thickness is about 160mm. 
End of article 
(i) Specify the real problem 


The article is concerned with the costs of heating a house and the possible 
benefits of insulation — but this does not state the problem sufficiently pre- 
cisely for the purposes of mathematical modelling. What is the actual problem 
that the article tackles? Use your own words. You should be able to give a 
satisfactory answer in a couple of sentences at most. 


(ii) Set up a model 


(a) There are many symbols in the passage: ¢,t,q, hy... Make a list of them 
all, explain precisely what each represents (‘t is thickness’ will not do!), 
give appropriate units of measurement, and say whether each represents 
a variable or a parameter. 


i] 


(6] 


(b) 


(d) 


The most important steps in the derivation of the formula for h(t) (Equa- 
tion 6) are those which lead to Equation 4 for the rate of loss of heat. 
Explain how the theory of heat conduction described in Unit 12 leads to 
Equations 2 and 3. Make sure that you state any assumptions about how 
the transfer of heat takes place that must be made to derive these equa- 
tions, whether or not they are mentioned in the article. Show in detail 
how Equation 4 follows from Equations 2 and 3. 


In the discussion of the cost of installation, the article states that it is 
easy to calculate the cost of the material because it will be proportional 
to the volume of material required. Show how this assumption leads to 
Equation 7 for i(t). 


Many assumptions needed to derive the model are mentioned in the ar- 
ticle, though they are not necessarily labelled as assumptions. Unfortu- 
nately, authors of articles on modelling sometimes overlook assumptions 
which in fact they use implicitly; and they sometimes state assumptions 
which are not really necessary. Bach of the following statements could 
possibly be an assumption which is needed for the development of the 
model. Decide, for each, whether or not it (or something equivalent) is 
explicitly mentioned in the article as an assumption in deriving the model; 
and decide whether or not it is strictly necessary to make that assumption 
in order to derive the model in the way described in the article. 


A The house has two stories ~ that is, the model would not apply to a 


bungalow. 
B All the heat lost into the roof space of the house while heating is 
required is supplied by the heating system in the house, 


C The insulating material used is fibye; 
D The rate of inflation is zero during the years over which the model 
is supposed to apply. 


ass. 


State one more assumption, mentioned in the article, which is necessary 
for the derivation of the model, and is not among those discussed in this 
part of the question (part (ii)), or in part (vi)(b) below. Use your own 
words; do not merely copy from the article, 


(iii) Formulate the mathematical problem 


(iv) 


The mathematical model obtained in the article is expressed in Equation 8, 
which gives the explicit formula for c(t) in terms of ¢. What would you regard 
as the main mathematical problem in this article? Is it 


Ori 


to calculate the yalue of ¢ using the data given in the final paragraph? 
to solve Equation 10 for t? 

to differentiate ¢(t) with respect to t? 

to find the value of ¢ for which ¢(t) is a minimum’? 


it something else entirely? 


Solue the mathematical problem 


(a) 


The right-hand sides of Equations 9 and 10 have been omitted deliberately 
from the version of the article printed above. Show in detail how the 
expression for t in Equation 11 is obtained from the formula for e(t) in 
Equation 8, and in the course of your answer give the missing parts of 
Equations 9 and 10. 


(b) It is stated after Equation 11 that this equation ‘gives the thickness of 


insulation for which the total annual cost is a minimum’, but nowhere is 
it shown that the stated value of t does correspond to a minimum (and 
not a maximum or a point of inflection), Show that the article is correct 
in asserting that it is a minimum. 


(3) 


(6) 


[2] 


(v) Interpret the solution 


(vi 


(a) One kind of qualitative check that can be made of the solution is to con- 
sider whether the formula for ¢ changes in the way one would expect if 
the values of the parameters are changed. This is something that can be 
investigated by considering the formula directly, or perhaps by sketching 
an appropriate graph. As an example, suppose we compare a house in 
the north of the UK with one in the south. We would expect it to be 
worthwhile to install more insulation in the north of the country, because 
winters are harder there than in the south. The model agrees with this ex- 
pectation, because a larger value of Aé@ leads to a larger value of ¢ for min- 
imum ¢(t). In each of the following cases, consider what happens to the 
value of ¢ for minimum e(t) according to the model, and comment on the 
extent to which this result corresponds to your expectations. You should 
assume in each case that all the other conditions remain unchanged. 


ty and oil goes up. 


A The cost of heating fuels such as gas, elect 
There is a rise in the interest rate. 
Heating is required for fewer days during the year. 
The price of insulating material comes down, 


mpoaw 


The insulating properties of the material available on the market for 
domestic loft insulation improve. 


(b) Does the model predict any circumstances in which it would make eco- 
nomic sense not to install insulation at all? Explain your answer. 


(c) With reference to the numerical calculation in the final paragraph, explain 
why 
A the U value of the ceiling material alone is 16 W m=? °C7!; 
B the cost of 1m* of fibreglass is £17.14; 
Cit is necessary to include the factor 1000 in the denominator of the 
fraction under the square root sign, 


Compare with reality 


It is not possible to compare this model with reality in the sense of comparing 
its predictions with what would happen in practice, without carrying out: some 
tests of the cost-effectiveness of insulation. But even so, there are some ques- 
tions that we can ask which help to evaluate how adequate the mathematical 
model is as a guide to installing loft insulation. 


(a) According to the article, the model predicts a layer of insulation thicker 
than is normally recommended. Does this, in your view, invalidate the 
model? Give reasons for your answer. 


(b) The outcome of the model clearly depends on the assumptions that are 
made in deriving it. So one important test is to consider how varying 
those assumptions would affect the results of the model. One assumption 
which may not always be realistic is that the insulation is laid by the 
householder. If this assumption is relaxed, it will be necessary to incor- 
porate labour costs in calculating i(t). Make the changes to the model 
which result from this, explaining clearly how you have taken account of 
the labour costs, and describe briefly what effects, if any, your changes 
will have on the outcome of the model. 


[2] 
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Question 2 


In this question you are asked to carry out some mathematical modelling for your- 
self. We shall first describe the situation, then make some suggestions about how 
you should tackle the problem, then give you some data, and finally explain what 
your tutor will be looking for in your solution, 


The situation 


One of the major motor-car manufacturers has been showing an intriguing TV 
advertisement for one of its models recently, This car is fitted with a driver’s 
airbag, and the purpose of the advertisement is to draw the viewer's attention to 
the benefits of having an airbag. In the advertisement, a car (of the relevant model) 
is driven very slowly off the flat roof of a tall building. It is seen falling more-or-less 
vertically downwards towards what appears to be a piece of white material on the 
ground where it would land. Just before the car is about to smash into the ground, 
however, the piece of material begins to move: it turns out to be a giant airbag, 
which inflates fully just in time to arrest the car, which lands on top of the bag at 
the exact moment when it is fully inflated, The bag then gently subsides, and the 
car is driven off without having suffered the slightest damage, 


Supposing that there was no electronic wizardry involved in the production of the 
advertisement, and that the events it shows actually happened as they appear to 
have happened, one cannot help wondering how the film-makers knew when to start 
blowing up the airbag, It would have been important to get the timing just right. 
If the inflation of the airbag had been left too late, then the pressure of gas in the 
bag wouldn’t have been sufficient to stop the car hitting the ground. On the other 
hand, the airbag begins to deflate automatically as soon as it has reached its fullest 
extent, so as to land the car gently (it behaves just like airbags inside cars in this 
respect, though they deflate in order not to obscure the windscreen any longer than 
is necessary). So if inflation had begun too early, the bag would also haye been too 
soft to arrest the car effectively. 


You are invited to use your mathematical modelling skills to solve the problem of 
when to start inflating the airbag. 


Method of attack 


You should develop a mathematical model of how the car falls which will allow you 
to predict how far above the ground it is at any given time after it has left. the top 
of the building. You will need to consider whether or not it is necessary to take 
any horizontal motion of the car into account, in addition to its vert motion. 
You will also have to decide whether or not to ignore air resistance, the possibility 
that the car will rotate as it falls, and any other similar factors. When you have 
modelled the falling motion of the car, you will have to work out how to use your 
model to predict when the car must hit the airbag if the ‘trick’ is to work properly. 


Data 

The following list should contain all the data you will need to solve the problem 

indeed, it should be more than enough. 

¢ The building is 100 metres tall. 

e The car is 4.8 metres long, 1.7 metres wide and 1.4 metres high; its mas: 
1300kg. 

= The acceleration due to gravity is 9.81ms7*. 

¢ Before it is inflated, the airbag covers an area of 1m?; when it is fully inflated, 
the area it covers is about 25m?. 

e =>. The time taken to inflate the airbag fully is 1 second. 


¢ = The airbag must have been inflating for at least 0.8 seconds for the pressure 
to have reached a level sufficient to support the car. 


e = The airbag will remain firm enough to support the car for a tenth of a second 
after it is fully inflated. 
e When the airbag is fully inflated, its top is 3 metres above the ground. 


Your solution 


You should write a brief account of your solution to this problem, Aim to keep to 
about the same length as the passage quoted in Question 1. That extract, and the 
drug therapy and the skid marks reports in Unit M, all provide useful examples of 
what is required. We suggest, however, that for definiteness you should base your 
solution as far as possible on the summary of Section 4 of Unit M, and that in your 
answer you should make sure you deal with the following points, 


(i) Specify the real problem 
Formulate in your own words a clear statement of the problem to be solved. {6] 
(ii) Set up a model 


Define precisely all of the parameters and variables you use, select symbols to 
represent them, and give appropriate units. (For example: ‘t is the time since 
the car left the roof, in seconds’.) Use symbols vather than numerical yalues 
for the quantities you need to formulate your model; this will make it easier 
for you to explain your model, and avoid mistakes. 6 


Carefully state the simplifying assumptions which you have made in formu- 
lating your model. (For example: ‘it is assumed that air resistance may be 
ignored.) 6 
Based on your assumptions, derive an expression for the downward component 
of the car’s acceleration while it is falling. Explain fully how you have obtained 
your result. 3 


(iii) Formulate the mathematical problem 


Explain how the knowledge of the downward component of the car's acceler- 
ation can be used to predict when it will be necessary to begin inflating the 
airbag. 4 
(iv) Solve the mathematical problem 


Obtain a formula giving the time at which the inflation of the airbag must 
begin, in terms of the height of the building, the height of the top of the 
airbag when it is fully inflated, and any other relevant parameters. 


(v) Interpret the solution 


Using the relevant data from the collection given above, obtain a quantitative 
answer to the initial question; that is, say how long after the car begins to fall 
inflation of the airbag should be started, 3 
Having obtained an answer, you must consider how reliable it is likely to be: 
if you were doing this for real, rather a large sum of money (the cost of a new 
car) would hang on the accuracy of your calculations. One obvious cause of 
possible error is the difficulty of measuring the height of the building sufficiently 
accurately: for example, the relative error in the height of the building could 
be 10%. What would be the effect on your answer of a 10% error in the height 
of the building? 4] 
Also consider briefly how sensitive your answer is to errors in the other data 
values. [3] 


(vi) Compare with reality 


We certainly do not suggest that you can or should carry out any experiments 
to test the reliability of your mathematical model. The car in the advertise- 
ment takes 14 seconds to fall from the top of the building to the inflated airbag. 
Do you think that your analysis throws any light on how the advertisement 
was really produced? Give reasons for your answer. [4] 


Identify one of the assumptions you made in part (ii) which you think could 
most. profitably be relaxed. Suggest how it could be improved, and say how 
you think that this will affect your answer. [4] 


EE 
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CMA MST204 44 Cut-off date 18 June 1998 


This assignment covers Units 15, 17, 18 and 19. 


Unit 15 


Questions 1 and 2 

A particle of mass 2kg moves with an acceleration, in ms~*, of 
2ti+t sin mtj, 

where i and j are Cartesian unit vectors, and ¢ is the time in seconds. 


1 What is the velocity of the particle, in ms~', at time ¢ = 1, if it is at vest at 
time t = 07 


Options 

Ai B 21 C 2i+nj D 
2 2. . 

E iris F Ut] G@ i+-j H 


2 What is the magnitude of the force, in newtons, acting on the particle at time 


t= ar 
Options 
A i B 2i Cc 4i D 4i+2j 
E 1 F 2 G4 H 2Vv5 


Questions 3 to 6 


A block of mass 2kg is in contact with a plane surface inclined at an angle of 30° 
to the horizontal. The coefficient of static friction between the block and the plane 
is equal to 0.59, whereas the coefficient of kinetic friction is equal to 0.56. 


Tn Questions 3 and 4 it is to be assumed that no forces act on the block other than 
its weight, the normal reaction and the frictional reaction. 


3 Decide whether or not the block can remain at rest on the inclined plane, 


Options 

A The block can remain at rest on the plane. 

B_ The block cannot remain at rest on the plane. 

4 — Select the option which is closest to the magnitude of the frictional reaction, 


|F;|, in newtons, if the block is moving down the plane, Here g denotes the 
magnitude of the gravitational acceleration, in ms~2. 


Options 
A 0.489 B 0.569 C 0.599 D 0.879 
E 0.979 F 1.009 G 1.029 H 1.739 


An additional force of magnitude P newtons is applied to the block in a horizontal 
direction when the block is at rest on the inclined plane, as shown in the diagram. 
This force is such that if it were any stronger, the block would begin to slide up the 
plane. 


5 Select the option which gives a correct expression for the magnitude of the 
normal reaction, |Fy|, in terms of P and g. 


6 Select the option which gives a correct expression for the magnitude of the 
frictional reaction, |Fy|, in terms of P and g. 


Options for Questions 5 and 6 
A g B V3g C P-g D -P+g 
E }P+V3g F -}P+V39 G iV3P-g H -1V3P +g 


Unit 17 


Questions 7 and 8 


Two blocks, of masses 5kg and 10kg, are connected by a light inextensible taut 
string. They move on a frictionless horizontal surface, along a straight line, A 
further light inextensible string is attached to the heavier block; it extends hori- 
zontally along the same straight line to a frictionless pulley, over which it passes. 
This second string is being pulled vertically downwards by a force of magnitude 
5 newtons, as shown in the diagram. 


Skg ~<— l0kg 


SN 


7 What is the magnitude of the common acceleration of the two blocks? 


Options 
A ims? B ims? Cc ims? D ims? 
E 2ms? F 3ms7* G 5ms™? H 10ms~? 


8 What is the magnitude T’ of the force exerted by the connecting string on 
either of the blocks? 


Options 
A ON B 14N D 2N 
E 2iN F 3N H 5N 


Questions 9 and 10 


A lorry is travelling along a straight horizontal road at a constant speed when a 
car, travelling in the same direction and at a constant but greater speed, comes up 
behind it and collides with it. As a result of the collision, the two vehicles become 
locked together. 


The mass of the lorry is 4000 kg and its speed before the collision is 15 ms~!, while 
the mass of the car is 1000 kg and its speed before the collision is 20ms7!. 

9 What is the speed of the centre of mass of the two vehicles before they collide? 
10 What is the common speed of the two vehicles immediately after the collision? 
Options for Questions 9 and 10 

A 15ms7! B 16ms~! C 16}ms7! D i7ms7! 


BE i7ims"! F 18ms"! G 19ms7! H 20ms* 


Questions 11 and 12 


Two particles collide at the origin. Before the collision, particle 1 is moving along 
the line y = # with speed u, while particle 2, whose mass is twice that of par 
is stationary. After the collision, particle 1 moves along the x-axis and particle 2 
moves along the y-axis, 


11 What is the speed of particle 2 after the collision? 


Options 

A 0 B hu Geis D 2u 

BE V2u F 2V2u iG Ee “ 

12 Is the collision elastic? 

Options 

A Yes B No 

Unit 18 

Questions 13 and 14 

Consider the fourth Taylor polynomial about « = —1 for the function e?", 


13 Which option gives the coefficient of (« + 1)* in the Taylor polynomial, correct 
to 3 significant figures? 

Options 

A 0.00564 B_ 0.0902 C 0.667 D 4.93 

E —0.00564 F —0.0902 G -0.667 H -4,93 


14 Select the option which is the error bound given by Taylor's Theorem for the 
fourth Taylor polynomial about x = —1 for e** in the interval —2 < x <0, 


Options 

A 3le B c 
D E F 
G H 


Questions 15 and 16 


An interpolating quadratic polynomial is used to approximate a function f(x) whose 
graph passes through the three points (—1,2), (0,1) and (1,3). 


15 Select the option which gives the interpolating polynomial. 


Options 

A l+a+2x? B 1-«+32* C 1-2¢+3: 
D 1+4u+ $2? E 1-$0+ $2? F -l+d2 
G -1+ het F207 H -4-Yr- $x? 


16 Suppose that |f@)(x)| <M in the interval [-1,1]. Select the option which 
is the best bound for the error function e(z) for the interpolating polynomial 
over the same interval. 


Options 

A |e(a)| < §la(@? -1)|M 
CC |e(a)| < ffa(a? - 1)|M 
E_ |e(x)| < $|x(a? —1)|M 
G e(x)| < d[2"|M 


le(w)| < 2yl(e — 1)(@ — 2)(@ — 3)|M 
le(w)| < §l(w = 1)(@ = 2)(@ - 3)|M 
e(a)| < $\(a — 1)(@ — 2)(@ — 3)|M 
le(w)| < p|a°|M 


mao 


Questions 17 and 18 


The approximate value of the definite integral 


2 
I cos(x") de 
Jo 


is evaluated numerically by the trapezoidal method and by Simpson's method using 
subintervals of width h = 0.5. 


17 Choose the option which gives the value, to four significant figures, obtained 
by the trapezoidal method. 

18 Choose the option which gives the value, to four significant figures, obtained 
by Simpson's method. 

Options for Questions 17 and 18 

A —0.7568 B 0.4615 C 0.4650 D 0.5271 

E 0.9300 F 1,832 G 2,790 H 3.665 


——— 


Unit 19 


Question 19 


The trapezoidal method is to be used, if possible, to obtain an approximate solution 
to the differential equation 


y! = (x+1)y—2? 
with the initial condition 
y(1) =1. 
Select the option which gives the approximate value of y(1.2), correct to 3 decimal 
places, if the step size is h = 0.2, or select Option G if you think that the trapezoidal 
method cannot be used. 
Options 
A 0,274 B 1.226 Cc 1.398 
D 1,622 E 1,723 F 1,851 


G_ The trapezoidal method, being an implicit method, cannot easily 
be used for this equation, 


Question 20 


If the step length h is reduced by a factor of 2 when using the Taylor series method of 
order 3, by what approximate factor is the local truncation error reduced, assuming 
that h is small? 


Options 
A 2 B 4 Cc 8 D 16 
E 32 F 64 G 128 H 256 


Question 21 


For what positive values of the step size h is the Taylor series method of order 2 
guaranteed to be absolutely stable for the differential equation 


f= = + cotx 
x 

with the initial condition 

y(1) = 2, 
over the interval 1 < a < 3? Choose the largest possible interval for h. 
Options 
A The method is unstable for all values of h. 
B h<02 Cc h<04 D h<06 
E h<l7 F h<20 G h<50 
H_ The method is stable for all values of h. 


Questions 22 and 23 
A numerical method uses the recurrence relation 
Your = ¥ + £A(3Y) + 2Y,,) 
to find approximate solutions for differential equations of the type 
y' = m(x,y). 
22 Select those TWO of the following statements which are TRUB. 
Options 
A The method is implicit. 
B_ The method is explicit. 
C_ The method is consistent with the differential equation. 
D The method is inconsistent with the differential equation. 
[There are TWO correct options] 


23 Select the option which is the principal term in the local truncation error for 
this method. 


Options 
A they! B 2h*y!! Cc snty! D Shey! 
B -dh?yt F =2h?y; G =py! H - Rhy 


—_—. 


Question 24 
The following methods can all be used to find numerical approximations for the 
solution of a certain first-order differential equation, with a given initial condition. 
The same step size is used in each case, and it may be assumed that with this step 
size, the method is absolutely stable. Which one of the methods would you expect 
to give the least accurate approximation? 
Options 
A Euler's method 

The Taylor series method of order 2 


B 
C_ The Taylor series method of order 3 
D_ The trapezoidal method 

E 


The predictor-corrector method (Euler-trapezoidal) 


aoe 
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This assignment covers Units 20, 21 and 22. 


Unit 20 


Questions 1 to 6 


Consider the matrices 


e-[o a], 9-[2 a], ®=[2 3): 
a=) 2 er x=["5]. y=(2 4]. 


1 Which one of the following options is equal to 2(P — R) +8Q? 


Options 
=8) 5 i 5 -9 5 
a | ‘| es cE + S [s | 
11 6 5 13 4 14 
D tio | a Ip ll E Le ‘| 


2 Which one of the following options is equal to Q-!? 


Options 
pou 1 
3 3 
ce ES Celia 
L-3 a 
1-2 i 3 no 
a Ei 2S [: 2] Fie 2A 


3 Which one of the following options is equal to Q7R? 
Options 


9 8 9 3 4 
at el) er 
5 4 5 2 8 
Dey eo | e (i 4 Bali i 


4 Which TWO of the following options represent matrix products which cannot 
be formed? 


so 


Options 
A PS B Qy G@ By” D xfPs 
E ys Foxy G yx H SQ 


[There are TWO correct options.] 
5 Which one of the following options is the value of det(PQ) ? 


Options 
A -6 Br =4 Gc —2 D oO 
E 2 F 2 G4 H 6 


6 There are two values of the number /& for which the matrix Q + $R+ kP is 
singular. Select the option which gives both values to three decimal places, or 
select Option H if you think none of the others is correct, 


Options 

A —0.667 and —0.250 B  —0.667 and 3.500 C —0.250 and 3.500 
D = —0,075 and 1,675 E 0,075 and 1,675 F 0,250 and 0.667 
G 0,250 and 3,500 H_ None of the other options is correct. 


Questions 7 to 9 


Consider the matrix 


and its inverse B= A~!. 


7 What is the value of det A? 


Options 
A= 8-2 ef Di 
Ba Pr 3 @4 H5 


8 Select the option which gives the value of the element bj of the matrix B. 
9 — Select the option which gives the value of the element boo of the matrix B. 
Options for Questions 8 and 9 

A —2 B -!I c -} DO 

E i 1 G2 


Questions 10 and 11 


New coordinate axes are defined in terms of given coordinate axes Oxy as follows: 
the new origin O' is at the point (—2,1) with respect to the axes Oxy, and the new 
axes O'x'y' are at a clockwise angle of 30° with respect to the axes Ory, 


10 Choose the option which gives the coordinates with respect to O'x'y! of the 
point whose coordinates with respect to Oxy are (—1,—3). 


11 Choose the option which gives the coordinates with respect to Oxy of the point 
whose coordinates with respect to O'x'y' are (1, —4). 


Options for Questions 11 and 12 


A (-2v3+4,4v3—4) 
Cee -ieey 

E (v5 -2 = a F 
(v3 


B 


i=) 


G 


w 
i 
& 
a 
+ 
aS 
+ 
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Question 12 
The matrices A, B and C are square and of the same size. and A is non-singular; 
Tis the unit matrix of the same size. Which TWO of the following statements are 
NOT necessarily true? 
Options 

If AB = 0, then B=0. 

If B’ = 0, then B=0. 


A 
B 

C If B? =0, then 1 +B is non-singular and (I+ B)~! = 
D Forall C, (1+C)? =1+20+C? 
E 

F 

G 

H 


For all B and C, (B + C)(B — C) = B? — C?. 

For all B, (B7)? = B. 

For all non-singular A, A” is non-singular. 

For all B and C, det(B7C”) = det(BC). 
[There are TWO correct options.] 


Unit 21 


Question 13 


A certain 2 x 2 matrix P has eigenvalues 
values of P? + 2P -1? 


—2 and 3. Which option gives the eigen- 


Options 
A —-9and 14 B -8 and 15 Cc -Tand 14 D -7and 16 
E -land 14 F -land 16 G and 15 H land 16 


Questions 14 to 17 


Consider the matrix 


rece 


14 The direct iteration method for finding eigenvalues is applied to this matrix, 
with initial vector yy =[0 1)". Which option gives the scaled vector y,? 


Options 
A [8 =77 B fs 7" e (=) =)" 
D [$ a)” E [-1 -7]” F ft #7’ 


T r 
Lee es 8 i 

e [ Vits Fl H [as Tl 

15 The inverse iteration method with p= 0 for finding eigenvalues is applied to 
this matrix with initial vector yp =(0 1)”. Which option gives the scaled 
vector y;” 

Options 

A [8 5)” B fl 


E [-4 5]? F fi 


16 The LR method for finding eigenvalues is applied to this matrix. Which option 
gives the matrix Lo? 


Options 
10 1 0 10 
«as ay ey 
41 81 41 
10 5 -8 5 0 
e lee] * bse ied 
17 Which option gives the matrix Up when the LR method is applied to the given 
matrix? 
Options 
5-8 5 =8 5 -8 
ata] = bo] ° bb Sl 
& =o 5 -8 5 -8 
> © f ae @ 


Question 18 


The matrix 


1 & 9.36: 
Aen Pee 8 
9 4 1 4 
To) 9! edt 


has eigenvalues approximately equal to —16.81, —3,32, —1.19 and 25.32. Which of 
the following options gives a value of p which could be used in the inverse iteration 
method in order to find a more accurate value for the eigenvalue near —3. 


Options 
A -10 B -2 GG D 2 E 10 


Unit 22 


Question 19 


“4 
‘The matrix eae has an eigenvalue \ = 1 with a corresponding eigenvector 
203 8 
{1 —1)", and an eigenvalue \ = 5 with a corresponding eigenvector [1 1)". 


Select the option which is the general solution of the system of differential equations 


& = 32, +220, 
dy = 2a, + 3x9. 
Options 
A B Cc ae! + be® 
—ae! + be 
D x, =ae! — be Ex, = 3ae! + 2he™* Fx; =ae** + be 
wv = ae’ + be aq = 2ae’ + 3be™ ay = ae” + be 


20 


Questions 20 and 21 


The matrix 
2 =2 
=2 32 
has an eigenvalue \ = 0 with corresponding eigenvector [1 1]”, and an eigenvalue 
\ =4 with corresponding eigenvector [1 —1]". 


20 Select the option which is the general solution of the system of differential 


equations 
iy = 2ary — 229, 
dq = —2a, + Qa. 
Options 
A a =ae Bx, =ae" Cc ne” + he 7 
ro = ae" £2 = —ae" ae” — be~** 
D E x, =at+be" F 2x, =at+be 
rt = —at + be" 
G H 


21 Select the option which is the general solution of the system of differential 
equations 
#, = 2x1 — 2x2, 
Hg = —2a, + 2x9, 


In the options, C,, C2, C3 and Cy are arbitrary constants. 


Options 

A |] =o,[!]+4e,]_?] 
wa | Pallas 2) =i] ~ 
xy 1 2 2 

B = ae + Cae") 
we =1 

ec {#]=- 1] (Cicos2e + Casin20) 
ing -1 

p |*l=a/?]4 1) (cae + Caen) 
ao) * {1 -1|? - 


5 
| 1 [1] + [1] (@cos2e + eusinay 


EO es he Jere | 1] (Cae + cue 
we 1 1 
ears 1] (ert ea + | 3] (Creos21 + crsin20 


Question 22 


The matrix 

-1 1 

-1 -1 
has an eigenvalue \ = —1+7% with a corresponding eigenvector [1 i)”, and an 
eigenvalue \ = —1 —i with a corresponding eigenvector [1 —]’’. 


Select the option which is the general solution of the system of differential equations 


&, = —2, +2, 
&: 2 — %. 
Options 
A x, = Ae‘ cost + Be'sint Bx, =Ae' cost + Be! sint 


xq = Ae' sint + Be' cost = Ae' sint — Be cost 


Cx, = Ae’ cost + Be‘ sint Dx, = Ae’ cost + Be'sint 
wg = —Ae' sint + Be’ cost Xq = —Ae' sint — Be' cost 
Ex, = Ae cost + Be“ sint Fx, = Ae cost + Be“ sint 
xg = Ae“ sint + Be~' cost wg = Ae~'sint — Be~' cost 
Ga = Ae cost + Be™ sint H x, = Ae™' cost + Be sint 
rg = —Ae~‘sint + Be‘ cost ay = —Ae' sint — Be™' cost 


Questions 23 and 24 


The matrix 


s( 


-2 -2 
=] -$8 


has an eigenvalue \ = —1 with corresponding eigenvector [2 —1]", and an eigen- 
value \ = —4 with corresponding eigenvector [1 1)”. A particle moves in a plane 
in such a way that its coordinates (w, y) satisfy 


#= Ar, 
where r=[x yJ?. 


23 Along which TWO straight lines can the particle move in simple harmonic 


motion? 
Options 
A Y=2z B y=2 


D y=-22 E y=-2" 
[There are TWO correct options.] 


22 


24 Select the option which gives the general form for the position of the particle, 
at time ¢, for any initial conditions. In the options, C;, C2, C3 and Cy are 
arbitrary constants. 


Options 


A «= 2C; cost + 2C2 cos 2t + Cz sint + Cy sin 2t 
y = —C; cost — C2 cos 2t + Cy sin t + Cy sin 2t 


B_ x= 2C, cost + 2C2sint + C3 cos4t + Cy sin 4t 
y = —C, cost — Cy sint + C3 cos 4t + Cy sin 4t 
C «= 2C\ cost + 2C2 sin t + Cy cos 2t + Cy sin 2t 
y = —C; cost — Cy sint + C3 cos 2t + Cy sin 2t 
D 2C;, cost + 2C2 sint + C3 cos 2t + Cy sin 2t 


ou 


—C\ cost — Cy sint — C3 cos 2t — Cy sin 2t 


eH 


E 2=Cycost + Cosint + 2C3 cos 2t + 2C, sin 2t 
y = C, cost + C2 sint — C3 cos 2t — C, sin 2t 


TMA MST204 05 Cut-off date 2 July 1998 


This assignment covers Units 17, 18, 19, 21 and 22. 


There are five questions in this assignment; Question 3 and Question 4 are both 
allotted fewer than 25 marks, however, and the total number of marks for the 
assignment is 100. 


Note that there is no tutor-marked assignment question on Unit 20; this unit is 
assessed by computer-marked assignment questions only. 


Question 1 (Unit 17) 
Note that the total number of marks available for this question is 25. 


A toy rocket is initially at rest on its launch pad and points vertically upwards, 
the motors being switched on at time t= 0. The total initial mass of the rocket 
(including its fuel) is 2kg. The fuel is burnt in such a way that the mass m of the 
rocket at time ¢ satisfies 

din 1 


pt Tet | «(OS tS 10). 


The rocket contains sufficient fuel for the motors to run for 10 seconds, and the 
exhaust gases are ejected vertically downwards at a constant speed of 490 ms~! 
relative to the rocket casing. For the motion of the rocket you may neglect any 
effects due to air resistance and assume that the magnitude g of the acceleration 
due to gravity is equal to 10ms~*. Throughout the question all the variables are 
measured in the appropriate SI units. 


(i) Find an expression for the mass m/(t) of the rocket for the period 0 <t< 10 
when its motors are running. 


(ii) Draw a diagram showing the external forces acting on the rocket during the 
period in which the rocket remains on the launch pad with the motors on, 


Find the upward reaction force exerted by the launch pad on the rocket at 
time ¢ during the period in which the rocket remains on the launch pad with 
the motors on. 


(iii 


Show that the rocket leaves the launch pad 2 seconds after the motors are 
switched on. 


(iv 


(v) Draw a diagram showing the external forces acting on the rocket during its 
vertical flight. 
(vi) Show that during its flight, the speed w of the rocket satisfies the differential 
equation 
a 980t 
dt — 200 —t? 
(vii) Find the speed of the rocket at the instant when the fuel is just exhausted, 
expressing your answer to an accuracy of one decimal place. 


-10  (2<t< 10). 


Question 2 (Unit 18) 
Note that the total number of marks available for this question is 25, 
(a) (i) Show that the quadratic polynomial p(x) which interpolates the function 
fle) =cosx 
at the points « = 0, 0.5 and 1 is 
p(x) = 1 — 0.029 97x — 0.429 732°, 


where the coefficients are correct to 5 decimal places, 


Sy 


(b) 


(ii) The equation 
cosa = 3x 


has a single root, which is near 0.3. By using p(x) as an approximation 
to cos, and solving a quadratic equation, find an approximation to this 
root better than 0.3, giving your answer to 3 decimal places. 


(iii) Use the Newton-Raphson method to determine the root of the equation 
cos x = 3x to an accuracy of 6 decimal places. 


Consider the differential equation 


dy _2_42 
ae ple = 
dx v (@ 20) 


with initial condition y(1) = 0. Use the Taylor series method of order 2 with 
step size h = 0,1 to find an approximation to y(1,2), where y is the solution 
to the equation with the given initial condition. Give your answer correct to 
6 decimal places. . 


Question 3 (Unit 19) 


Note that the total number of marks available for this question is 13. 


This question is concerned, once again, with the solution of the differential equation 


dy 
ae 


-y* 


with the initial condition y(1) = 0. 


(i) 


(ii) 


(iii) 


Starting from the general form of the recurrence relations for the predictor- 
corrector (Buler-trapezoidal) method, write down the recurrence relations 
for finding a numerical solution to the above differential equation using this 
method, Use your calculator to find an approximation to y(1.2), where y is 
the solution satisfying the given initial condition, using a step size h = 0.1. 
Give your answer correct to 6 decimal places. 


‘The values shown in the table below were obtained for the approximations to 
y(2) given by the predictor-corrector (Buler-trapezoidal) method with various 
step lengths h by using a computer. Plot a graph of these approximations to 
y(2) against h?. 


Step length h Approximation to y(2) 


0,05 1.565 785 
01 1.564497 
0.125 1.563 497 
0.2 1,559 028 
0.25 1.554916 


What does your graph indicate to be the relationship between the predictor- 
corrector (Euler-trapezoidal) approximations to (2) and the step length h for 
small h? Use your graph to estimate the true value of y(2) correct to four 
decimal places. 


By using one of the entries in the above table and the estimate for the true 
value of y(2) you obtained in part (ii), calculate the maximum value of the 
step length i which you would have to use in the predictor-corrector (Buler- 
trapezoidal) method to ensure an accuracy of six decimal places in the approx- 
imation to y(2). 


(4) 


(4 


w 
& 


Question 4 (Unit 21) 
Note that the total number of marks available for this question is 12. 


Find the eigenvalues and corresponding eigenvectors of the matrix 


1 4-1 
-1 0 -1]. 1) 
—2 =4 0 


Question 5 (Unit 22) 
Note that the total number of marks available for this question is 25, 
This question is concerned with the system of differential equations 
a = Tay — 3x2 —€**, 
ta = 62, — 2x9. 
(i) Express the system of differential equations in the form 
X(t) = Bx(t) + h(t), 


where 
x(t) = fer(t) a(t], h(t) = [Pa(t) a(t)” 
and B is a constant 2 x 2 matrix. [2] 
(ii) Find the eigenvalues and corresponding eigenvectors of the matrix B. Hence 
find a matrix P such that P~'BP is a diagonal matrix. [8] 
(iii) Find the general solution of the system of differential equations. {15] 


Printed in the United Kingdom by 
‘The Open University 


26 


